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g i ABSTRACT 
OO 

We report on an observation of the low-mass X-ray binary 4U1626-67 
^ | performed during the BeppoSAX Science Verification Phase. An absorption 

feature at ~ 37 keV, attributable to electron cyclotron resonance, has been 
discovered in its pulse averaged spectrum. The inferred neutron star magnetic 
field strength is 3.2 • (1 + z) x 10 12 G, where z is the gravitational redshift. The 
feature is deep and narrow and is resolved in both the broad-band fit and in the 
ratio of observed counts to those seen from the Crab. The cyclotron resonance 
energy is in good agreement with the empirical relation between cyclotron 
energy and high energy cutoff, while its width is in agreement with the expected 
Doppler broadening of thermal electrons at the cyclotron resonance frequency. 
The broad-band 0.1-200 keV spectrum is well fit by a two-component model: a 
0.27 ± 0.02 keV blackbody and a power law with a photon index of 0.89 ± 0.02. 
This is the first broad-band observation made after the change from spin-up 
to spin-down that occurred in mid 1990: it confirms the harder spectrum with 
respect to those observed in the 2-10 keV range. 

Subject headings: binaries: close — stars: neutron — pulsars: individual: 
4U1626-67 — X-rays: stars 
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1. Introduction 



Following the first detection in Her X-l flTrumper et al. 1978| ), line features in the 
10-110 keV energy range were detected in some X-ray pulsars ( jMakishima fc Mihara 
1992|; |Grove et al. 1995|) . These features are interpreted as being due to electron 
cyclotron transitions in the ^ 10 12 G magnetic field of the neutron star (jMakishima 



|fc Mihara 19921) . Because the energy of the fundamental cyclotron harmonic, E cyc , is 



related to the neutron star magnetic field strength in units of 10 12 G, B\2, by the relation 
E cyc = 11.6 B12 ■ (1 + z)^ 1 keV, where z is the gravitational redshift, the observation of 
cyclotron resonance features (CRFs) in pulsar spectra gives a direct measurement of the 
magnetic field of the neutron star. From a theoretical point of view CRFs are expected to 
be visible as absorption features. This was shown e.g. by Nagel (1981) for any reasonable 
optical depth. The feature is not strictly due to absorption, but rather to scattering of 
photons in the wings of the line. Up to now CRFs have been observed mainly in the spectra 
of young high-mass binary systems, with the notable exception of Her X-l flTrumper et al. 
19781) . 



The 7.7 s X-ray pulsar 4U1626-67 is one of the few pulsators among LMXRBs. Its 
timing history is characterized by a sudden change in its spin state that occurred on 1990 
June QUhakrabarty et al. 1997| ). Before this date the source was steadily spinning-up, while 
from the transition to the present time 4U1626-67 is spinning-down. Besides coherent 



pulsation, 4U1626-67 also shows quasi-periodic oscillations (QPOs) both in X-rays ( (Shinoda 
|et al. 1990| ) and in the optical QChakrabarty 199~8| ), with a centroid frequency of ~ 40 mHz. 



The pulse-phase averaged spectrum of 4U1626-67 has been described in terms of a 
two component model: an absorbed blackbody with a temperature kT ~ 0.6 keV and a 
power law (|Kii et al. 1986| ). The spectrum also exhibits a high energy cutoff at ~ 20 keV, 
interpreted as due to a possible CRF at that energy ( |Pravdo et al. 1979|) , and an emission 
Ne line complex at ~ 1 keV ( |Angelini et al. 1995[ ). While the spin transition did not affect 
the high energy part of the spectrum, within the accuracy of previous measurements, it did 
affect the 2-10 keV part ([Vaughan fc Kitamoto 19"97|) : the power law photon index changed 
from ~ 1.6 ( Pravdo et al. 19"79| ) to <^ 0.7 QAngelini et al. 1995|) , with a correspondingly 
significantly lower 2-10 keV flux, that shows an overall fading by a factor four from the 
HEA01/A2 measurement ( |Chakrabarty et al. 1997] ). 
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2. Observations 

The BeppoSAX satellite is a program of the Italian Space Agency (ASI), with 
participation of the Netherlands Agency for Aerospace Programs (NIVR), devoted to X-ray 
astronomical observations in the broad 0.1-300 keV energy band. The payload includes 
four Narrow Field Instruments (NFIs) and two Wide Field Cameras. The NFIs consist of 
Concentrators Spectrometers (C/S) with 3 units (MECS) operating in the 1-10 keV energy 
band ( [Boella et al. 19971 ) and 1 unit (LECS) operating in 0.1-10 keV flParmar et al. 1997|) , 



a High Pressure Gas Scintillation Proportional Counter (HPGSPC) operating in the 3-120 
keV energy band ( [Manzo et al. 1997] ) and a Phoswich Detection System (PDS) with four 
scintillation detection units operating in the 15-300 keV energy band ( Frontera et al. 1997] ). 



During the Science Verification Phase a series of well known X-ray sources were 
observed in order to check the capabilities and performances of the instruments on-board 
BeppoSAX . 4U1626-67 is one of these sources and it was observed by BeppoSAX on 1996 
August 6, from 14:41 to 21:50 UTC, and from 1996 August 9 00:10 to August 11 00:00. 
Data were telemetred in direct modes, which provide complete information on time, energy 
and, if available, position for each photon. The net exposure times for the four NFIs were 
31, 97, 43 and 52 ks for LECS, each of the MECS units, HPGSPC and PDS, respectively. 
These differences are due to the LECS only being operated in satellite night-time, rocking 
of collimators (for HPGSPC and PDS) and different filtering criteria during passages in the 
South Atlantic Geomagnetic Anomaly (SAGA) and before and after Earth occultations. 
Furthermore, HPGSPC data from the August 6 observation were not included because of 
telemetry problems. The spurious spikes (At ^ 2 s) that affect the PDS data were also 
filtered. These spikes occur at a frequency of few per orbit, can reach several hundreds 
counts per second and are likely due to the local environmental particle background. Their 
spectrum is generally soft, showing a strong cutoff at ~ 40 keV (Guainazzi fc Mattcuzzi 



1997] ) . Their presence can severely affect both a spectral and timing analysis, expecially for 



moderately faint sources in the PDS bandpass, such as 4U1626-67. 

The data analysis of all the NFIs but the LECS was carried out in the XAS V2.0 
framework flChiappetti fc Dal Fiume 19971) . For the LECS data the SAXLEDAS Vl.4.0 
was used ( [Lammers 1997] ) . We used the 1997 September revisions of the response matrices 



for all NFIs but LECS, for which we used the matrix generated by LEMAT V3.5.3. 
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3. Spectral Analysis 



The analysis of the LECS data has been already presented flOwens, Poster broek, fc 



[Farmar 1997| ), and will not be discussed here: we will include their result on the Ne line 
complex. For the imaging instruments, the standard extraction radius of 8' was used to 
derive the source events. The rocking collimated instruments (HPGSPC and PDS) net 
counts were obtained with the standard background subtraction procedure, described in 
details in Manzo et al. (1997) and Frontera et al. (1997). The systematic error in PDS 
background subtraction is consistent with zero for observations longer than 100 ks and 
negligible for sources whose flux is higher than 1 mCrab, such as 4U1626-67 flGuainazzi fe 
Matteuzzi 1991 ). 



3.1. The X— ray continuum 

The spectra of each BeppoSAX NFIs were first fit separately in order to derive the 
best model able to describe the broad-band spectrum. We find that a model including low 
energy absorption, a blackbody, a power law and a high energy cutoff adequately describes 
all the NFIs spectra but the PDS, yielding xl s of 1.13 (81 degrees of freedom, dof) for 
LECS, 1.11 (540) for MECS, 1.29 (144) for HPGSPC, and 1.89 (56) for PDS. Most of 
the contribution to the high xl f° r PDS is due to deviations at ~ 35 keV. We used this 
model to describe the broad-band 0.1-200 keV spectrum. We included in the fit a variable 
normalization for each of the NFIs (we treated separately each MECS unit), in order to 
take into account for the known calibration uncertainties between the instruments. As high 
energy cutoff we used the standard X-ray pulsar model exp[(£' c — E)/Ef], where E c and Ef 
are the cutoff and folding energy, respectively flWhite, Swank, fc Holt 1983| ). The smoother 
Fermi-Dirac cutoff ( [l'anaka 1986| ) did not adequately describe the high energy tail. After 
the inclusion of the Ne line complex at 1 keV we obtained a best fit with a reduced xl °f 
1.369 for 706 dof. The fit results are summarized in Table [I]. The normalizations among 
NFIs (MECS2 as reference) are 1.04 for LECS, 1.05 for HPGSPC and 0.72 for PDS. 

Our data do not show the presence of an Iron K-shell line in 6.4-6.9 keV: the 3a 
upper limit on its equivalent width is 21 eV, slightly more stringent than the 33 eV value 
obtained by ASCA ( [Angelini et al. 1995| ). The total 0.1-200 keV X-ray luminosity is 
7.7 x 10 34 erg s" 1 d£ pc . The fluxes in the bands 2-10, 10-60, and 10-200 keV are 1.7, 4.4, 
and 4.5 x 10~ 10 ergs cm -2 s _1 , respectively. 
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3.2. Cyclotron resonance 

From the analysis of the residuals we were led to add an extra component to our 
spectral model. Both a Lorenzian flMihara 1995[ ) or a Gaussian in absorption ( [Soong et al. 



1990|) at ~ 35 keV, interpreted as a CRF, improved the fit, yielding xl s of 1.146 and 1.232 
for 703 dof, respectively. An F-test shows that the improvement is significant at 99.99%. 
Adding this component improved the fit to the 15-200 keV PDS spectrum, yielding a xl of 
1.256 (53 dof; probability of chance improvement 2.4 • 10 -5 ). The Gaussian line FWHM 
of 7 ± 2 keV is consistent with the PDS energy resolution at the CRF energy, ~ 17.5% 
( |Frontera et al. 1997| ). Figure |1] shows the best-fitting spectral model with the residual CRF 
Gaussian line profile obtained by setting the line normalization to zero. The residuals for 
each NFIs are also shown. 

The two CRF models are equivalent: the probability of chance improvement from the 
Lorenzian to the Gaussian model is 34%. The different Ecrf found with the two models 
is an effect due to the modeling: the Lorenzian model approximates the fall-off of the 
spectrum by increasing the cutoff energy, and shifting the cyclotron energy to values lower 
than those obtained by an absorption Gaussian (compare, in Table [j], the best fit continuum 
parameters with the CRF line included with those obtained from the continuum without 
the line). This is a known effect, already observed for the CRF in Her X-l ( pal Fiume 
et al. 1998Q and Vela X-l ( |Orlandini et al. 1998| ). 



Following the suggestion by Pravdo et al. (1978) of a possible CRF at ~ 19 keV, we 
tried to include in our broad-band model a fundamental CRF at this energy (and therefore 
to consider the ~ 35 keV CRF as the second harmonic). This attempt was unsuccessful, 
both with and without the high energy cutoff component. 

In order to better characterize the 4U1626-67 CRF we examined the ratio between the 
PDS and Crab count rate spectra. This ratio has the advantage of minimizing the effects 
due to the detector response, and the uncertainties in the calibrations. As we can see from 
the upper panel in Fig. 0, there is an evident change of slope at ~ 35 keV. The same ratios 
performed on the Her X-l ( Pal Fiume et al. 1998|) and Vela X-l ( Prlandini et al. 1998|) 



spectra show a change of slope at ~ 40 and ~ 57 keV, their well known cyclotron energies. 

To enhance the effect due to the CRF, we multiplied the 4U1626-67/Crab Ratio by 
the functional form of the Crab spectrum, i.e. a simple power law with index equal to 
2.1. In this way we emphasize deviations of the 4U1626-67 spectrum from its continuum 
without making any assumption on its form. The result is shown in the second panel of 
Fig. where it is evident that deviations occur only at E ^ 30 keV. 

Finally, in the lower panel we show the ratio between the previous function and the 
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4U1626-67 continuum model, with a, E c and Ef taken from Table |l|. In this way we 
enlarge all the effects due to line features, although we introduce a model dependence. The 
absorption feature at ~ 38 keV is evident. This feature can be fit between 25 and 50 keV 
with both a Gaussian (shown in the figure) and a Lorenzian. The results are shown in 
Table ^| and are in agreement with those obtained from the broad-band spectral fit shown in 
Table [l]. The difference in the ocrf values is the effect due to the intrinsic energy resolution 
of the PDS instrument. The inferred magnetic field strength at the neutron star surface is 
therefore (3.2 ± 0.1) x 10 12 • (1 + z) G, where z is the gravitational redshift. 



4. Discussion 

This is the first 4U1626-67 broad-band spectrum obtained after the spin transition. 
From our flux measurement we can see that the overall bolometric flux of the source 
decreased by a factor 4.1 since the HEA01/A2 measurements. We find a power law index 
slightly higher than those found in narrow-band spectra ( |Angelini et al. 1995| ; Pwens, 



Poster broek, fc Parmar 1997| ). This is an effect due to the broad-band fit, in which we 



simultaneously take into account both the low-energy power law and the high-energy 
cutoff. 

Our 4U1626-67 CRF measurement, together with its cutoff energy, fits the empirical 
relation found in X-ray binary pulsars between these parameters: higher the cyclotron 
resonance energy, higher the cutoff — and therefore harder the spectrum ( [Makishima et al 



|1990| ; [Frontera fc Dal Fiume 1989| ). Our CRF measurement also is in agreement with the 



expected electron Doppler broadening. Indeed, at the cyclotron resonance frequency ou c , 
electrons at rest absorb photons of energy huj c . For moving, thermal electrons the Doppler 
broadening Aud is predicted to be flMeszaros 1992| ) 

(2kT e \ l '\ „, 
Auj d = uj c \ - cost* (1) 

\m e c 2 J 

where kT e is the electron temperature, and m e c 2 is the electron rest mass. The angle 6 
measures the direction of the magnetic field with respect to the line of sight. Outside 
the range uj c ± Aojd the cyclotron absorption coefficient decays exponentially, and other 
radiative processes become important. From Eq. [l] and the CRF parameters given in 
Table [I] we obtain a lower limit to the electron temperature of ~ 5 keV, in reasonable 
agreement with the calculations of self-emitting atmospheres of Harding et al. (1984) , 
according to which the temperature for a column or slab with optical depth ~ 50 g cm" 2 is 
4-8xl0 7 °K, corresponding to kT e ~ 3.5-7 keV. 
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Finally, we want to mention that if we assume the QPO frequency as due to the 
beating between the pulse frequency and the Keplerian motion at the magnetospheric 
radius, then the 4U1626-67 magnetic field strength is related to its luminosity by the 
relation B\i ~ 5.5V-^37 ( [Shinoda et al. 199C ), where L 37 is the X-ray luminosity in units 
of 10 37 erg s _1 . Assuming a source distance 5 ^ d kpc ^ 13 flChakrabarty 1998|) , this 
corresponds to 2.4 <^ Bu 6.3, in good agreement with our measurement. 



The authors wish to thank the BeppoSAX Scientific Data Center staff, the Nuova- 
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Fig. 1.— The broad-band 0.1-200 keV spectrum of 4U1626-67 observed by the BeppoSAX 
NFIs. Left panel: count rate spectrum (crosses) and best-fitting continuum model 
(histogram) with the Ne complex and absorption Gaussian components (parameters from 
Table 1). The Gaussian normalization has been set to zero and the residuals from the model 
are shown. Right panel: detailed NFIs residuals from the broad-band best-fit model (for 
clarity only one MECS unit is shown). 

Fig. 2. — Upper panel: Ratio between the PDS 4U1626-67 and Crab count rate spectra. 
Note the change of slope at ~ 35 keV. Middle panel: the 4U1626-67/Crab ratio multiplied 
by E -21 , the functional form of the Crab spectrum. Deviations from the continuum (dotted 
line) occur for E ^ 30 keV. Lower panel: ratio between the (4U1626-67/Crab ratio x E -2 ' 1 ) 
and the model of the 4U1626-67 continuum as derived from the broad-band fit. A Gaussian 
fit to the CRF (dotted line) is also shown. 



Table 1. Best-fit spectral parameters a 



Parameter Value 







No Line 


Gaussian 


Lorenzian 


N H 


(10 21 cm' 2 ) 


1.1 ±0.2 


1.2 ±0.2 


1.1 ±0.2 


kT 


(keV) 


0.28 ±0.01 


0.29 ±0.01 


0.30 ±0.01 


r-bb 


(xd kpc km) 


1 c;+0.3 


1.5 ±0.2 


1.4 ±0.2 


a 




0.87 ±0.01 


0.86 ±0.01 


0.83 ±0.01 


E c 


(keV) 


21.1 ±0.4 


19.6 ±0.5 


29 ±2 


E/ 


(keV) 


7.8 ±0.3 


10.6 ±0.7 


9± 2 


i pow 




1.06 ±0.02 


1.05 ±0.02 


1.01 ±0.02 


ENe 


(keV) 


1.05 c 


1.05 c 


1.05 c 




(keV) 


0.04 c 


0.04 c 


0.04 a 


EW Ne 


(keV) 


48 c 


48 c 


48 c 


EcRF 


(keV) 




37 ±1 


33 ± 1 


CCRF 


(keV) 




3± 1 




EWcrf 


(keV) 




15 ±3 




TCRF 








1.6 ±0.3 






1.369 (706) 


1.232 (703) 


1.146 (703) 



a Uncertainties at 90% confidence level for a single parameter. 

b Flux is in units of 10~ 2 photons cm~ 2 s _1 at 1 keV. 

c The Ne line complex parameters are taken from Owens et al. 
(1997). 



Table 2. Fit results on the 4U1 626-67 CRF a 



CRF Parameter Gaussian Lorenzian 



E CRF (keV) 38.6 ± 0.9 38.1 ± 0.9 

<r crf (keV) 5.0 ±0.7 8 ± 1 

xl (dof) 1.048 (13) 2.725 (13) 



a Uncertainties are given at 90% 
confidence level. 
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